Single fluxon in double stacked Josephson junctions: Analytic solution. 
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We derive an approximate analytic solution for a single 
fluxon in a double stacked Josephson junctions (SJJ's) for 
arbitrary junction parameters and coupling strengths. It is 
shown that the fluxon in a double SJJ's can be character- 
ized by two components, with different Swihart velocities and 
Josephson penetration depths. Using the perturbation the- 
ory we find the second order correction to the solution and 
analyze its accuracy. Comparison with direct numerical sim- 
ulations shows a quantitative agreement between exact and 
approximate analytic solutions. It is shown that due to the 
presence of two components, the fluxon in SJJ's may have an 
unusual shape with an inverted magnetic field in the second 
junction when the velocity of the fluxon is approaching the 
lower Swihart velocity. 

Magnetic field, H , parallel to layers, penetrates stacked 
Josephson junctions (SJJ's) in the form of fluxons. Flux- 
ons in SJJ's are different both from Abricosov vortices in 
type-II superconductors, since they don't have a normal 
core, and from Josephson vortices in a single Josephson 
junction (JJ's), since magnetic field is spread over many 
JJ's. Despite of a large body of works concerning vor- 
tex related properties of layered superconductors, for re- 
view see e.g. Ref. 0, exact solution for a single fluxon in 
SJJ's is still missing. Distribution of magnetic induction, 
B, at large distances from a fluxon was studied in Refs. 
Pj and (|], for tightly packed stacks with identical and 
nonidentical layers, respectively. In Ref. Q a solution 
for weekly coupled double SJJ's was obtained using per- 
turbation approach. However, this solution is not valid 
for strongly coupled SJJ's, which is the most interesting 
case. The electromagnetic coupling is strong when inter- 
layer space periodicity of the stack, s, is (much) less than 
the effective London penetration depth, A a {,, as it is the 
case for high-T c superconductors. 

Recently, a simple approximate analytic solution for 
double SJJ's was suggested j5j. The approximate solu- 
tion was shown to be in a good agreement with numer- 
ical simulations in the static case for arbitrary coupling 
strengths and parameters of the stack. However, the ac- 
curacy of the solution was not strictly proven. Another 
question which is still opened is how the shape of the 
fluxon is changed in the dynamic case. This is a crucial 
question since fluxon dynamics determines perpendicular 
(c-axis) transport properties of SJJ's. 

In this paper we rigorously derive the approximate ana- 
lytic solution for a single fluxon in a double SJJ's, which 
is valid for arbitrary junction parameters and coupling 



strengths. Using the perturbation theory we find the 
second order correction to the solution and analyze its 
accuracy. Comparison with direct numerical simulations 
shows a quantitative agreement between exact and ap- 
proximate analytic solutions. It is shown that the fluxon 
in a double SJJ's can be characterized by two compo- 
nents, with different Swihart velocities and Josephson 
penetration depths. We also studied the transformation 
of the fluxon shape with increasing propagation veloci- 
ties. It is shown that due to the presence of two com- 
ponents, the fluxon in SJJ's in the dynamic case may 
have unusual shape with an inverted magnetic field in 
the second junction at high propagation velocities. 

We consider a double SJJ's with the overlap geometry, 
consisting of JJ's I and 2 with the following parameters: 
J c i -the critical current density, ti- the thickness of the 
tunnel barrier between the layers, di and Xsi - the thick- 
ness and London penetration depth of superconducting 
layers. Hereafter the subscript i on a quantity represents 
its number. The strength of electromagnetic coupling of 
junctions is determined by the coupling parameter, S, 
which varies from to 0.5 for the stack with identical 
JJ's. More details about definitions can be found in Ref. 
II . We consider frictionless fluxon motion with a con- 
stant velocity, u. The fluxon will be placed in junction 
1. 

The physical properties of SJJ's are described by cou- 
pled sine-Gordon equation (CSGE) |6|-|| for gauge in- 
variant phase differences, tfi. The problem with solving 
CSGE is the coupling of nonlinear sin ((ft) terms. To de- 
couple the variables, we take linear combination of equa- 
tion in CSGE and rewrite it as: 

^,2^-«n(^i,2) = -Bn,2(0, (!) 
where £ = x — ut is the self-coordinate of the fluxon, 

-Fl,2 = fl ~ «2,l¥>2, (2) 



Erifl = sin(tpi) - n 2 ,isin(ip 2 ) 

~ -K 2 ,1<P2 (1 - COs(lfi)) - 



- sin(ipi - «2,i<£2) 
O (^) • (3) 



Here \\ 2 = Af 2 (l — u 2 /c^ 2 ) and coefficients Kx$, 
characteristic Josephson penetration depths, A12, and 
characteristic Swihart velocities, c\ %2 are given by 
Eqs. (17,20,21) from Ref. [||. Such choice for coefficients 
Ki t 2 minimizes Er\,2 far from the fluxon center. The 
phase differences should satisfy boundary conditions: 
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^i(-oo) = 0,99i(0) = 7r, (px(+oo) = 2tt; 
^(±00,0) = 0. (4) 

Therefore, functions Er\$ have a form of ripple around 
zero value and will be considered as perturbation. Solu- 
tions of Eq.(l) can be easily found. Solutions of uniform 
Eq.(l), i.e. with zero r.h.s., are: 



Aarctan 



exp 



(?/Ai, 2 ) 



(5) 



From Eq.(2) we obtain the first approximation for 
phase differences: 



¥>2 = 



KlFl - K 2 F 2 
Kl - K 2 

Fi — F 2 

Kl — K 2 ' 



(6 a) 
(6 b) 



which coincides with the approximate analytic solu- 
tion, obtained in Ref. j5j. The approximate solution is 
asymptotically correct at large distances from the fluxon 
center and has correct values at x = 0, as follows from 
Eqs. (3,4). 

Next, we look for a solution of nonuniform Eq.(l) in 
the form Fi t 2 = Fi,2 + SFi^ 2 , where SFi y2 are corrections 
due to perturbation terms Er±2- 



cosh 2 



(e/Ai, 2 ) 



SF h2 = En, 2 , (7) 



with boundary conditions 8Fi y2 (±oo, 0) = 0. The so- 
lution of Eq. (7) is 



8 Fx , 2 = 01,2(0/1,2 + h,2(09i,2, 



(8) 



where 

/l,2 



l/cosfc(£/Ai,2), 



ffi,2 = sink — — 

\Al,2 / 



Ai,2Cos/i(^/Ai, 2 ) 
are partial solutions of the uniform Eq.(7) and 
-1 r« 



(9) 



ai 2 = 



2Ai,2 Jo 



bi, 2 



1 



2Ai, 2 



Er 1<2 (x')gx,2 {x')dx', 



Eri, 2 {x')fi, 2 (x')dx'. 



(10) 



The perturbation corrections to the approximate 
fluxon solution, Eq.(6) are 



Sifi = 
Sip 2 = 



K1SF1 - k 2 SF 2 

Kl - K 2 

8F1 - 8F 2 

Kl - K 2 



(11a) 
(lib) 



The corrections 5ipi, 2 are used to improve the approx- 
imate solution Eq.(6) and to estimate it's accuracy. 

Lets observe that from Eq. (3), KiEn ~ k 2 Et 2 . 
From Eqs. (9-11) it can be seen that for A2 = Ai, 
K1SF1 — k 2 SF 2 and 5ip\ = for arbitrary tp 2 ^ 0. More- 
over, taking into account Eq.(6 b), it can be shown that 

In Fig. 1, central part of a fluxon is shown for a 
stack of two identical JJ's with strong coupling, S'=0.495, 
for the static case, u—Q. Phase distribution in the full 
scale is shown in Fig. 3. Parameters of the stack are: 
dis = h,2 = O.OIAji, A S i_ 3 = O.IAji, where X n is 
the Josephson penetration depth of the single junction 1. 
Solid lines represent results of direct numerical integra- 
tion of CSGE, Eq. (1), dashed lines show the approx- 
imate analytic solution, Eq.(6), and dotted line shows 
the corrected analytic solution, tp 2a + S(p 2a , Eqs. (6, 11), 
in junction 2. Solid and dashed curves, marked as 5ipi, 2 , 
represent the overall discrepancy between numerical and 
approximate analytic solutions and the perturbation cor- 
rection, Eq.(ll), respectively. 

From Fig. 1 it is seen that correction to the fluxon 
image, ip 2 , in the second junction vanishes far from and 
in the fluxon center, while the accuracy decreases at dis- 
tances ~ Aji from the center. Such behavior is expected 
from the shape of perturbation functions Eri y2 in Eq. 
(1). On the other hand, for junction 1, correction 5ip\ 
is small in the whole space region and analytic solution 
gives an excellent fit to the "exact" numerical solution, in 
agreement with discussion above. The most crucial test 
for the approximate solution is the accuracy of derivative 
at x — 



Vi(0) _ «i6i(0)/Ai-«a6a(0)/A a 



Vi(0) Ki/Ai-Ka/Aa 
An estimation for u = yields 

2aKiK 2 (Xi ~ A 2 ) 2 (\i 1 + A^ 1 + ^&-) 



(12) 



Vi(0) 



(13) 



where a is a factor of the order of unity and Ao given 
by Eq.(26) from Rcf. Q. From Eq. (13) it is seen that 
both Sf'i (0)/<p[ (0) and it's derivative with respect to A2 
goes to zero at A2 = Ai, as discussed above. 

In Fig. 2, the maximum of Sipi (x) (top panel) and the 
relative correction to derivative at x = 0, Sip'i (0)/ip'i(0), 
(bottom panel) for u=0 are shown as a function of 
Jc2/Jci for four different coupling parameters S*=0.495 
(solid lines), ,5=0.433 (dashed lines), 5=0.312 (dotted 
lines) and 5=0.127 (dashed-dotted lines). For 5=0.495 
parameters of the stack are the same as in Fig.l; for 
5=0.433 d t = 0.5X Si ; for 5=0.312 d t = X Sl ; and for 



2 



5=0.127 di = 2Xsi- From Fig. 2 it is seen that the ac- 
curacy of solution improves with decreasing S. This is 
naturally explained by a decrease of (f2, see Eq.(3) and 
decrease of splitting between Ai^, see Eq. (20) from Ref. 
||. However, even for strongly coupled case, the analytic 
solution, Eq.(6), gives quantitatively good approximation 
not only for the value, but also for the derivative of ipi 
for arbitrary parameters of the stack. The gray solid line 
in the bottom panel of Fig. 2 shows the estimation for 
Vi(0)M(0)> calculated from Eq.(13) for 5=0.495. It 
is seen, that estimation gives qualitatively correct result. 
Namely, 6<pi (0)/<pi (0) vanishes both for J C 2/Jci — > 0, as 
Jc2/Jci, and for J c ij Jc\ — ► oo, as yj J c i/ J C 2- Note, that 
for J c 2/J c i 0, Sipi vanishes even though the splitting 
of Ai,2 becomes extremely large ||. 

Sofar we have considered the static case, u—0. On 
the other hand, according to Eq. (6), radical changes 
should take place in the dynamic state, and the shape 
of the fluxon in SJJ's may become qualitatively differ- 
ent from that in the single Josephson junction. Indeed, 
as the velocity approaches the lower characteristic ve- 
locity, u — ► Ci, Ai — > 0, i.e. the F\ component of the 
fluxon contracts, while contraction of the F2 component 
remain marginal. This implies, that at u — ► c\ , the fluxon 
in SJJ's consists of a contracted core and uncontracted 
"tails" decaying at distances much larger than the core 
size. Such behavior is clearly different from that in a sin- 
gle Josephson junction. We note that characteristic ve- 
locities, Ci ; 2 may depend on u, therefore, contraction of 
each component, Fi,2, may be different from Lorentz con- 
traction. Another interesting consequence of the approx- 
imate analytic solution, Eq.(6), is that with increasing u, 
the magnetic field in the second junction may change the 
sign with respect to that injunction 1. Such behavior was 
predicted in Ref. from the approximate analytic solu- 
tion, Eq. (6), and it was suggested, that this will result 
in attractive fluxon interaction in SJJ's. The magnetic 
induction in SJJ's is equal to H 



B 1 = 



Bo 




2(l-S 2 ) 



(14 a) 
(14 b) 



where Hq = 7rA $0 A and A 1,2 are defined in Ref. M. 

From Eq.(14) it is seen that estimation of magnetic 
induction in SJJ's requires the accuracy of derivatives 
in both junctions, while Eq.(6) is not valid with the ac- 
curacy of derivative at x = for if2- Therefore, more 
elaborate analysis is needed for the study of magnetic 
field distributions in the fluxon. 

In Fig. 3, we show a) phase distributions and b) mag- 
netic field distributions in the fluxon for different fluxon 
velocities. Parameters of the stack are the same as in 
Fig. 1. Solid lines in Fig. 3 a) represent results of di- 



rect numerical simulations of CSGE, Eq.(l), and dotted 
lines show the approximate analytic solution, Eq.(6). It 
is seen that quantitative agreement between " exact" and 
approximate solutions sustain up to c\. For the case of 
identical junctions, considered here, exactly one half of 
the fluxon belongs to each of the components, F\^- In- 
deed, from Fig. 3 a) it is seen that for u ~ c\ there 
is a contracted core at x = with a one n step in ipi. 
On both sides of the core, there are two tt/2 tails, which 
are slowly decaying at distances ~ A2 3> Ai. Solid and 
dashed curves in Fig. 3 b) represent numerically simu- 
lated profiles, B lt 2 (x), injunctions 1 and 2, respectively. 
From Fig. 3 b) it is clearly seen, that with increasing 
fluxon velocity, a dip in B2 develops in the center of the 
fluxon. At velocities close to ci, B2 (0) changes sign and 
finally at u = Ci, B2 (0) = —B\ (0). Such behavior is 
in agreement with predictions of Ref. || . From numeri- 
cal simulations we have found, that the fluxon shape in 
double SJJ's is well described by Eq.(6) up to at least 
u w 0.98ci for any reasonable parameters of the stack, 
although the accuracy of the approximate solution may 
decrease with increasing u. The decrease of accuracy is 
caused by the increase of if 2 as is seen from Fig. 3. In 
this case, perturbation correction, Eq. (11), should be 
taken into account. 

In conclusion, a simple approximate analytic solu- 
tion for a single fluxon in a double stacked Josephson 
junctions for arbitrary junction parameters and coupling 
strengths is derived. It is shown that the fluxon in a 
double SJJs can be characterized by two components, 
with different Swihart velocities and Josephson penetra- 
tion depths. Using the perturbation theory we find the 
second order correction to the solution and analyze its 
accuracy. Comparison with direct numerical simulations 
shows a quantitative agreement between exact and ap- 
proximate analytic solutions for all studied parameters 
of the stack and fluxon velocities up to at least 0.98ci. 
It is shown that due to the presence of two components, 
the fluxon in SJJ's may have an unusual shape with an 
inverted magnetic field in the second junction at large 
propagation velocities. This may lead to attractive fluxon 
interaction in the dynamic state of SJJ's 
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FIG. 1. Central part of a fiuxon in a stack of two identical 
JJ's with strong coupling, 5=0.495, for the static case, u=0. 
Solid and dashed lines represent "exact" numerical and ap- 
proximate analytic solutions, respectively. Dotted line shows 
the corrected analytic solution, ip 2a + S(f2 a - 

FIG. 2. Perturbation correction to the approximate ana- 
lytic solution is shown as a function of Jc'ij Jc\ for the max- 
imum of (5y>i (x) (top panel) and the relative correction to 
derivative at x = 0, (V'i(0)/(p'i(0), (bottom panel) for u=0 
and for four different coupling parameters. The gray solid 
line in the bottom panel shows the estimation from Eq.(13) 
for 5=0.495. It is seen that the analytic solution, Eq.(6), gives 
quantitatively good approximation not only for the value, but 
also for the derivative of ipi for arbitrary parameters of the 
stack. 

FIG. 3. Fiuxon shapes in double SJJ's for different fiuxon 
velocities, w/ci=0, 0.61, 0.92, 0.98, 0.998, 0.9999 (from left 
to right curve). In Fig. 1 a) phase differences, <£>i,2, are 
shown. Solid and dotted lines represent "exact" numerical 
and approximate analytic solutions respectively. Fig.l b) 
shows magnetic inductions Bi,2 obtained numerically. The 
existence of contracted and uncontracted components and the 
sign inversion of £2(0) at u ~ ci is clearly seen. 
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